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We analyze the breakup of invariant tori in Hamiltonian 
systems with two degrees of freedom using a combination of 
KAM theory and renormalization group techniques. We con- 
sider a class of Hamiltonians quadratic in the action variables 
that is invariant under the chosen KAM transformations, fol- 
lowing the approach of Thirring. The numerical implemen- 
tation of the transformation shows that the KAM iteration 
converges up to the critical coupling at which the torus breaks 
up. By combining this iteration with a renormalization, con- 
sisting of a shift of resonances and rescalings of momentum 
and energy, we obtain a much more efficient method that 
allows to determine the critical coupling with high accuracy. 
This transformation is based on the physical mechanism of the 
breakup of invariant tori. We show that the critical surface 
of the transformation is the stable manifold of codimension 
one of a nontrivial fixed point, and we discuss its universality 
properties. 



I. INTRODUCTION 

The existence of invariant tori plays a significant role 
for the long-time stability of Hamiltonian systems. The 
KAM theorem JJ^] states that the tori with frequency 
vectors which satisfy diophantine conditions are sta- 
ble under small perturbations. Conversely, it has been 
shown H that for large perturbations these tori do no 
longer exist. For two-dimensional systems, Greene 
formulated a criterion that allows to determine the exis- 
tence of a KAM torus by analyzing the properties of a 
sequence of nearby periodic orbits, namely the resonances 
whose winding ratios are the continued fraction approxi- 
mants of the winding ratio of the considered torus. When 
the amplitude of the perturbation is at its critical value, 
these resonances open gaps in the torus and it breaks up 
into a Cantor set (Aubry-Mather set) |p[-|lT|. 
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In order to study the self-similar scaling 
served for the breakup of invariant tori |lz 
ization group ideas were proposed for two-dimensional 
area- preserving maps |^|,[h|. For Hamiltonian systems 
with 1.5 degrees of freedom, Escande and Doveil [ p7|p8| 
set up an approximate renormalization scheme that com- 
bines KAM transformations with a rescaling of phase 
space. 

The idea of renormalization group analysis for Hamilto- 
nian systems is to construct a transformation TZ as a gen- 
eralized canonical change of coordinates acting on some 



space of Hamiltonians such that the iteration of TZ con- 
verges to a fixed point. If the perturbation is smaller 
than the critical one, TZ must converge to some Hamilto- 
nian Hq for which the equations of motion show trivially 
the existence of a torus of a given frequency vector. All 
Hamiltonians attracted by this trivial fixed point have 
an invariant torus of that frequency: this statement can 
be considered as an alternative version of the KAM the- 
orem Jl9| . If the perturbation is larger than critical, the 
system does not have a KAM torus of the considered fre- 
quency and the iteration of TZ diverges. The domain of 
convergence to Hq and the domain of divergence are sep- 
arated by a surface invariant under the action of 1Z. The 
main hypothesis of the renormalization group approach 
is that there should be another nontrivial fixed point on 
this critical surface that is attractive for Hamiltonians 
on that surface. Its existence has strong implications 
concerning universal properties in the mechanism of the 
breakup of invariant tori. The analysis of the renormal- 
ization for area-preserving maps fl5| gives support to the 
validity of this general picture. The aim of the present 
work is to give similar support to this picture for Hamil- 
tonian flows. The main ideas had been announced in Ref. 

The transformation we define (KAM-RG) has two main 
parts: a KAM iteration which is a change of coordinates 
that reduces the size of the perturbation from e to e 2 , 
and a renormalization transformation which is a com- 
bination of a shift of the resonances and a rescaling of 
momentum and energy. It acts within a space of Hamil- 
tonian systems with two degrees of freedom, quadratic in 
the action variables. An essential aspect of the present 
approach, based on a formulation of the KAM theorem 
by Thirring plf , is that the KAM and renormalization 
transformations we use map this space into itself. In or- 
der to analyze the strong coupling regime and to reach 
the critical coupling, the KAM-RG transformation has to 
converge at least up to the critical surface at which the 
torus breaks up. In fact, we show that the KAM iteration 
as well as the KAM-RG transformation converge all the 
way to the critical coupling. Numerically, the KAM-RG 
transformation is a much more efficient method to de- 
termine the critical coupling. The analysis of the KAM- 
RG transformation shows that the critical surface is the 
stable manifold of codimension one of a nontrivial fixed 
point. 

We construct the KAM transformation by two alterna- 
tive methods: by the Lie transformations and by trans- 
formations defined by a generating function. The moti- 
vation to use two different transformations is two-fold. 
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First we verify that both approaches lead to the same 
results, and the Lie transformation is more efficient for 
numerical implementation. Further, although the two 
transformations lead to quantitatively different nontriv- 
ial fixed points, they have the same critical exponents, 
in accord with the general ideas of the renormalization 
group approach. 

In Sec. It we describe the renormalization transforma- 



(1.5) 



where C L = {v E 1? \\v x \< L, \v 2 \ <L}. We define (/) 
the mean value of / by 



(/> = 



T 2 (2tt) 2 



(1.6) 



tion. In Sees. |llj and [IV], we define the KAM iteration where T 2 = [0, 2tt] x [0, 2tt]. In the following sections, we 

df 

will use the notation df= — — for any function of the 

dip 

angles. 



of the transformation by the two methods. In Sec. [v|, we 
give our numerical results, and in particular, we show ev- 
idence of the existence of an even nontrivial fixed point. 
In Sec. VI, we describe the behavior of the KAM-RG 



transformation when odd perturbations are included. 

We consider the following class of Hamiltonians with 
two degrees of freedom, quadratic in the action variables 
A = (Ai,A2) and described by three scalar functions of 
the angles ip = (ipi,ip 2 )- 

H(A, ip) = ^m(<p)(Sl • A) 2 + [«„ + g(<p)Sl] ■ A + f(ip) , 

(1.1) 

where Wo is the frequency vector of the considered torus, 
and n = (1, a) is some other constant vector, not parallel 
to uiq. This class of Hamiltonians has been considered 
by Thirring pl| in its non-degenerate version 

H(A, ip) = • M(<p)A + [u + g(<p)] ■ A + f(<p) , 

(1.2) 

where M is a 2 x 2 matri x su ch that detM ^ 0, and g a 
vector. The Hamiltonian (1.1) is such that 



det 



d 2 H 
dAdA 



0. 



(1.3) 



i.e. it does not satisfy the twist condition, but the KAM 
theorem is also valid under this condition (se e p^ |). The 
advantage of the family of Hamiltonians l \l . 1| ) in the 
present context is that they are characterized by three 
scalar functions of the angles and a con stant a, instead 
of six functions for Hamiltonians (1.2). This allows a 
more precise numerical treatment of the problem. The 
essential features are already contained in the space of 
Hamiltonians (1.1). In particular, the nontrivial fixed 
point one obtains starting with (1.2) is of the form (1.1). 
The functions to, g, f are represented by their Fourier 
series e.g. 



(1.4) 



The numerical implementation of the transformation re- 
quires a truncation of the Fourier series. We will approx- 
imate / by 



II. RENORMALIZATION TRANSFORMATION 

We construct the KAM-RG transformation combin- 
ing two parts: a KAM transformation (to, <?, /, a) > 
(to', g' , /' , a) and a renormalization consisting of a shift 
of the resonances and a rescaling of the actions and of 
time (m',g',f',a) <— > (to", g", /", a'). The renormaliza- 
tion scheme described in this section is for a torus of 
frequency u> = (1/7,-1) where 7 = (1 + v5)/2, but 
this scheme can be adapted to quadratic irrationals. 
The KAM-RG transformation is composed of four steps: 
1) a KAM transformation which is a change of coordi- 
nates that eliminates terms of order 0(e), where e is the 
size of the perturbation; this transformation produces 
terms of the order 0(e 2 ) and does not change $7 = (1, a) 
(see Sees. 



Ill and 



IV) 



2) a shift of the resonances: a canonical change of coordi- 
nates that maps the next pair of daughter resonances of 
the sequence of rational approximants into the two main 
resonances, 

3) a rescaling of energy (or equivalently of time), 

4) a rescaling of the action variables (which is a general- 
ized canonical transformation). 

The aim of this transformation is to treat one scale at the 
time. The steps 2), 3) and 4) are implemented as follows: 
The two main resonances (1,0) and (1,1) are replaced 
by the next pair of daughter resonances (2, 1) and (3, 2), 
i.e. we require that cos[(2, 1) • ip'] = cos[(l,0) • <p"] and 
cos[(3, 2) • ip'] — cos[(l, 1) • ip"]. This change is done via 
a canonical transformation (A',ip') (N~ 2 A', N 2 ip') 
with 



N 



2 1 
1 1 



This linear transformation multiplies u>q by 7 -2 (since Wo 
is an eigenvector of N); therefore we rescale the energy 
by a factor 7 2 in order to keep the frequency fixed at u)q. 
A consequence of the shift of the resonances is that £1 is 
changed into ft' = (1, a'), where a' = (a + l)/(a + 2). 
Then we perform a rescaling of the action variables: we 
change the Hamiltonian H' into 



l^.ip 



Phys. Rev. E 57, 1536 (1998). 



3 



with A such that the mean value (m") is equal to 
1. Since the rescaling of energy and the shift N 2 
transform the quadratic term of the Hamiltonian into 
7 2 (2 + a) 2 m'((f')(n' ■ A') 2 /2, this condition leads to 
A = 7 2 (2 + a) 2 {m!). This condition has the following 
geometric interpretation in terms of self-similarity of the 
resonances close to the invariant torus: The rescaling 
magnifies the size of the daughter resonances, and places 
them approximately at the location of the original main 
resonances. This can be seen by the following heuristic 
argument: In order to estimate the position of the reso- 
nances, we assume that e is small and that m' w (m'). 
The equations of motion for H' are 



A' ps const, 

ip' ps (m')(n ■ A')fi + u:o. 



(2.1) 
(2.2) 



The position of the resonance v is given by the condition 
v ■ tp' = 0, i.e. it is located at A' such that 



U -A' re — 



1 LJq, ■ V 

m') n ■ v 



(2.3) 



The linear change of coordinates TV 2 gives the new posi- 
tion of the resonance v 



SI' -A' 



1 



Wo ■ v 



(m')7 2 (2 + a) 2 fl' ■ v 



(2.4) 



Thus, the rescaling of the actions A' i— > A" = A'/A with 
A = 7 2 (2 + a) 2 (m') places the resonances at the location 
of the original resonances. 

In summary, the renormalization rescales m, g, f and 
O = (1, a) into 



m"(cp) 



m' (N~ 2 tp) 



g"( V )= 1 2 (2 + a)g> (j\rV) , 
f"(<p)= 1 *(2 + a) 2 (m , )f (J\T V) , 
, 1 + a 



Q 



(2.5) 

(2.6) 
(2.7) 

(2.8) 



The iteration of the transformation (2.8) converges to 
a* = 7 _1 . It means that Q converges under successive 
iterations to f2* = (1,1/7), which is orthogonal to ujo 
and is the unstable eigenvector of iV 2 with the largest 
eigenvalue 7 2 . 

We remark that this renormalization scheme can also 
be implemented on a more general class of Hamilto- 
nians quadratic in the action variables considered by 
Thirring |l|: 



H(A,tp) = ±A.M(<p)A- 



[wo +»(¥>)] A 



f(f) , 
(2.9) 



where M is a 2 x 2 matrix and g a vector. This class 
of Hamiltonians is also invariant under KAM transfor- 



mations (see Sees. [II and IV). The renormalization de- 
scribed above changes the direction of g and M. The 



vector g is renormalized into 71(g) = j 2 N 2 g. The it- 
eration of 1Z converges to the unstable eigenvector of 
N 2 : g — > g£l*. The matrix M is renormalized into 
1Z(M) = N 2 MN 2 /{N 2 MN 2 ) n . This transformation 
has only one stable fixed point M* = 17* ® $7*. Thus 
the iteration of the renor maliz ation transformation on 
Hamiltonians of the form ( |2.9| ) (which can satisfy the 
twist condition or not) converges to a twistless Hamilto- 
nian (1.1) with f2 = $"2*. 



III. LIE APPROACH TO THE KAM 
TRANSFORMATION 

The Poisson bracket of two functions of cp and A is 
given by 

r , , df dg df dg 
{L9}= d^-dA-dA-W (3 - 1} 

We will work with Lie transformations Us '■ {f, A) 1— » 
(V',A') 

A' = e^^'^A = A-{S,A} + {5, A}} . . . (3.2) 



f-{S, V } + - [ {S,{S, V }} 



(3.3) 



generated by functions S linear in the action variables, 
of the form 



S{A, ip) = Y{tp)n ■ A + Z(ip) + aQ ■ , 



(3.4) 



characterized by two scalar functions Y , Z , and a con- 
stant a. The expression of the Hamiltonian in these new 
variables is obtained by the following equation 



H'(Ai, V ') = e + ^H(A,<p)\ (A ,, v , ) 

= H + {S,H} + ^{S,{S,H}} 



(3.5) 



A consequence of the linearity of S in A is that the Hamil- 
tonian H' is again quadratic in the actions, and of the 
form 



ff'(AV) = ^V)(tt-A') 2 



[wo + g'(<p')Sl] ■ A' + f'(tf') 



(3.6) 



We notice that ft is unchanged by this transformation. 
Following the approach of Thirring pl| , we consider the 
scalar functions g and / of order 0[s), and m of or- 
der one. We determine S such that g' and /' are of 
order 0(e 2 ). The nth iteration of this transformation 
will produce g and / of order 0(e 2 ) and m of order 
one. The idea is that m(ip) does not need to be elim- 
inated. In order to show the existence of a torus of 
frequency u>q that is located at A — 0, it suffices that 



the iteration reduces the Hamiltonian (1.1) into one with 
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/ = 0, g = but m((p) ^ 0. This is an immedi- 
ate consequence of the equations of motion associated 
to H = m(ip)(fl ■ A) 2 /2 + u - A. The fact that m{<p) 
does not need to be eliminated is what allows to work 
with canonical transformations that are linear in the ac- 
tions. This has the important practical advantage that 
the KAM transformation leaves invariant the sub space of 
Hamiltonians quadratic in the actions [Eq. (1.1) or Eq. 

©]■ 

The expressions of <?' and /' up to the order 0(e 2 ) are 

g'( l p')=g(<p')+u; -dY 

+m{<p')(n-dZ + an 2 ) + 0{e 2 ), (3.7) 
/V) = f(«p') +u, -dZ + a^-n + 0{e 2 ). (3.8) 



Thus in order to eliminate the terms of order 0(e), we 
determine: 

1) Z(tp) such that the term independent of the action 
variables is of the order 0(e 2 ). The function Z must 
satisfy the equation 



/ + u>o ■ dZ = const, 
which has the solution 



(3.9) 



(3.10) 



The mean value of Z is not determined by Eq. (3.9). We 
choose it equal to zero. 

2) Next, we determine a and Y(tp) such that the linear 
term in the action variables becomes of the form [u>o + 
0(e 2 )fl} ■ A'. This leads to the condition 



g + iv ■ dY + mtl ■ dZ + man 2 = 0, 

which has the solution 

(g) + (mSl ■ dZ) 



and 
Y(<p) 



E 



LJ ■ V 



fl 2 (m) 



(g u + (mfl ■ dZ) v + m„af2 2 



(3.11) 



(3.12) 



(3.13) 



The transformed Hamiltonian (3.6) is constructed by 
defining ffM = H and for i = 1, 2, ... by the re- 
cursive relation 

H« +1 \A,<p) = {S(A,<p),H®(A,<p)} 



+g^(<p)n-A + f^ +1 >(<p), (3.14) 



( V )(n-Af 



(»+!)( 



which leads to 



00 



(3.15) 



This can be expressed in terms of the image of the three 
scalar functions (m,g,f) given by the following equa- 
tions: 



v i=0 i=0 i=0 



(m<°W 0) ,/ (0) ) = (rn, <?,/), 

m {1) = 2mO • dY - Yfl ■ dm, 
g {1) = gfldY-Yfl-dg 



(3.16) 
(3.17) 

(3.18) 



+mfl ■ dZ + mafl 2 + w ■ dY, (3.19) 
= -Yfl ■ df + gU ■ dZ + gaVL 2 + u; • dZ, (3.20) 

TO (i+1) = 2m W fi dY-YVt- dm {€} , (3.21) 
g (i+i) = g (*)n.BY-YH-dg® 

+m w O- dZ + m {l) an 2 , (3.22) 

f {l+1 ^ = -YU ■ df [i) + g {l) Vt ■ dZ + g {i) afl 2 , (3.23) 

for i > 1. 



IV. GENERATING FUNCTION APPROACH TO 
THE KAM TRANSFORMATION 

In this section we describe another construction of the 
KAM transformation. The transformation is now taken 
as a canonical transformation Up : f tp, A) i— > (tp' , A') 
defined by a generating function p5| characterized by 
two scalar functions Y , Z, of the angles, and a constant 
a, of the form 

F(A\ V ) = {A' + an)-cp + Y(cp)n ■ A' + Z(cp), (4.1) 

leading to 

dF 

A = — = A 1 + n ■ A'dY + atl + dZ, (4.2) 
dF 

<p' = — = <p + Y(<p)n. 



dA> 



(4.3) 



Inserting Eq. (4.2) into Eq. (1.1), we obtain the expres- 
sion of the Hamiltonian in the mixed representation of 
new action variables and old angle variables 



H(A',^)= 1 -m^){n-A') 2 



>o + g(<p)n]-A' + f(<p) , 



with 



m = (1 + ft ■ dYfm, 

g = g + lj ■ dY + mb + ft ■ dY (g + mb) , 
f = f + LJ -dZ+ X -mb 2 + gb, 



(4.4) 



(4.5) 
(4.6) 

(4.7) 
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where b((p) = aft 2 + fl • dZ. We notice that the KAM 
transformation does not change fl. 



We determine the generating function (4.1) such that TLo 
Up vanishes to the first order in e. This leads to the 
conditions 



wo ■ 9Z + f 



const, 
m (aft 2 



n ■ dz) = o. 



(4.8) 
(4.9) 



We recall that the functions g and / are of order 0(e) 
and to is of order one; as a consequence Y , Z and a are of 
order 0(e). We notice that these equations are the same 
as Eqs. ( |3.9[ ) and ( |3.11 ) which determine the generator of 
the Lie transformation. The present tr ansf ormation and 
the Lie transformation described in Sec. Ill are canonical 
transformations with identical linear part [i.e. O(e) part] 
but different nonlinear terms (of the higher order in e). 
The main practic al d ifference with the Lie transforma- 
tion is that Eq. (fD|) which de term ines the new angles 
has to be inverted. Equations (4.1 



and ( |4.9| ) are solved 
by representing them in Fourier space. They define the 
generating function F as 



Z{ V ) 



f v e 



IV-if 



(g) + (mft ■ dZ) 



fl 2 {m 
i 



Y( V ) = — Gfc + (mil • BZ\ 



-m v aft 2 



Thus the scalar functions of H become 

to = (1 + ft ■ dY) 2 m, 
g = -u) Q ■ dYCt ■ dY, 

J=\{g-uo-dY) ( a n 2 + n-dz) 



(4.10) 
(4.11) 

(4.12) 



(4.13) 
(4.14) 

(4.15) 



We notice that ifi, g and / are given by products and 
sums of functions whose Fourier coefficients are explic- 
itly known. We expand these functions in Fourier series 
e.g. fh(tp) = J2v 



The expression of the Hamiltonian in the new angle vari- 
ables requires the inversion of Eq. fl4.3| ) . The Jacobian of 
this transformation is 



dct 



dipi 



= |i + n-0Y| 



(4.16) 



The scalar functions ttt,', g' and /' are respectively to, g 
and / expressed in the new angle variables e.g. m!(ip') = 
m[(p((p')]. The Fourier coefficients of to' are determined 
by the following integrals 

/„;, = l ^-rn'^e-^', (4.17) 



With the change of variables ip' i— > ip, we can write 



(Pip 

T2 (2^)* 



det 



d(pi 



fh{<p)e- iu <v +Y ^ n \ 

(4.18) 



Therefore m' v , g' v and f v can be expressed in the following 
way: 

m v — y ' fh u iC v i v , g v = S ' g v iC v i u , 

fl = zZf» ,c ^> ( 4 - 19 ) 

I/' 

where 

Cy, v = \ —^\l + n-dY\e l( "'-"> v e- lv - nYM . 

(4.20) 

In order to compute C v i v , we choose the Gauss quadra- 
ture that approximates an integral as a sum over a lattice 
of constant step: 



T 2 (2tt) z 



.... J— y 

— ( M+1 ) 2 , i =t'. > M 



lim 



G 



2n 



2tt 



M + \ M + 1 



(4.21) 



This quadrature is exact if all the non-zero Fourier modes 
of the considered function G are inside Cm- For instance 
in the case of the computation of C„v for the identity 
transformation (Y = 0), since v' — u G C2L, one needs 
to take M > 2L. For the general case, if we expand 
the function to be integrated as a power series of e, we 
notice that the algorithm with M = (k + \)L gives an 
approximation up to 0(e k ). For instance this algorithm 
with k = 5 gives an accurate approximation of the exact 
transformation and allows to compute all its properties. 



V. DETERMINATION OF THE CRITICAL 
COUPLING; FIXED POINT OF THE KAM-RG 
TRANSFORMATION 

We start with the same initial Hamiltonian as in Refs. 

H(A,cp) = i(n- A) 2 +uj q - A + ef(<p) , (5.1) 

where £1 = (1,0), w = (1/7,-1), 7 = (1 + V^)/2, and 
a perturbation 



f(tp) = cos(i/i • ip) + cos(z^2 • ip), 



(5.2) 



where ui = (1,0) and V2 — (1,1). We perform an 
iteration of the KAM transformation described in Sec. 
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[II or in Sec. IV. The two methods give qualitatively 



the same results. The algorithm using the Lie transfor- 
mation is numerically more efficient. The following re- 
sults are those obtained by the Lie transformation. Wc 
represent all the functions by their Fourier series trun- 
cated by retaining only the coefficients in the square Cl 
which contains (2L + l) 2 Fourier coefficients. For fixed 
L we take successively larger couplings e and determine 
whether the KAM iteration converges to a Hamiltonian 
with / = 0, g — 0, or whether it diverges (/, g — > oo). By 
a bisection procedure, wc determine the critical coupling 
e c (L). We repeat the calculation with larger numbers of 
Fourier coefficients, to obtain a more accurate approxi- 
mation. In Fig. 1, we show s c {L) i.e. the dependence 
of the critical coupling on the number of Fourier coeffi- 
cients retained. We observe that e c (L) decreases with L 
in a stepwise manner. It stays essentially constant except 
at the values of L where a new rational approximant of 
the frequency u>o is included, corresponding to a reso- 
nance at the next smaller scale. The size of the jumps 
diminishes approximately geometrically, and we can ex- 
trapolate to obtain the value e c (L) — > 0.0276. This value 
is close to the critical coupling e c — 0.0275856 obtained 
by the Greene criterion |5],^6| which is surmised to yield 
the exact value. This gives a numerical evidence that the 
KAM iteration can be expected to converge in the whole 
domain of existence of the torus. 

In Figs. 1 and 2, we show the values of the critical cou- 
pling e c (L), calculated by the KAM-RG transformation 
which is a combination of a Lie transformation (Sec. 
Ill) and a renormalization transformation (Sec. ||). We 
obtain e c £ [0.027585,0.027595], which is in very good 
agreement with the value e c = 0.0275856 obtained with 
the Greene criterion. We observe that the KAM-RG 
transformation gives very high precision already with few 
Fourier coefficients e.g. e c (L = 5) = 0.0276633. 
The improvement with respect to the KAM iteration is 
not only quantitative; the disappearance of the steps is 
a strong evidence that the KAM-RG transformation we 
have constructed captures the essential physical mecha- 
nism of the breakup of the tori. 

By iterating the KAM-RG transformation starting from 
a point on the critical surface, we observe that the pro- 
cess converges to a nontrivial critical point 7J* , which we 
characterize by the Fourier coefficients of the three func- 
tions /*,<?*, m* and f2* = (L7" 1 ). Figure 3 shows the 
weight of the Fourier coefficients of /* . We observe that 
the nonzero coefficients are strongly concentrated on a 
band around the direction D± perpendicular to the line 
of resonances. D± is the expansive direction of the map 
v i— ► N~ 2 v, i.e. the direction of the frequency vector Wq. 
The decrease of the size of the coefficients along D± is 
quite slow. The Fourier coefficients of g* and m* have a 
similar overall behavior, but they decay faster in the D± 
direction (see Figs. 4 and 5). By linearizing the KAM-RG 
transformation around the fixed point H*, we calculate 
the critical exponents. There is only one with modulus 
greater than one, denoted by 5. This implies that the 



critical surface, which is the stable manifold of is of 
codimension one. The value we obtain for the relevant 
critical exponent is 6 £ [2.67, 2.68] which is quite close to 
the one obtained by MacKay for area-preserving maps 
]l5t (S — 2.65), and to the one obtained by Escandc et 
at. with the approximate scheme (6 — 2.75) p7|| . 
For the scaling factor at the nontrivial fixed point, we 
obtain numerically A„ = 17.944, which is very close to 7 6 
(the scaling factor at the trivial fixed point). This value 
can be compared with the one given for area-preserving 
maps A* = 18.827 obtained in Refs. 0||||]§ . 
We remark that if we consider the more general starting 
Hamiltonian of the form ([T|) H = A-MA/2 + iv Q -A + f 
with the twist condition detM ^ 0, the KAM-RG itera- 
tion has the same nontrivia l fix ed point as we have found 
for twistless Hamiltonians (1.1). 

We have also performed this analysis using the KAM 
transformation constructed with a generating function 
as described in Sec. |v|. The results (using fewer Fourier 
coefficients than the Lie transformation) are qualitatively 
similar. The nontrivial fixed point is quantitatively dif- 
ferent, but the critical exponents are the same. This is 
what is usually expected in a renormalization group ap- 
proach. 



VI. SYMMETRIES OF THE 
TRANSFORMATION: GENERALIZATION TO 
NON-EVEN PERTURBATIONS 

In the previous section, we found numerically the exis- 
tence of a nontrivial fixed point for the KAM-RG trans- 
formation starting with an even initial perturbation on 
the critical surface 



/o(v) = -Mcos(u! ■ ip) - Pcos{v 2 ■ <p), 



(6.1) 



where M, P > and u i, i>2 correspond to the two main 
resonances. For a perturbation containing also odd terms 
sin(i/i • cp) and sin(i/2 • tp), the KAM-RG transformation, 
acting on the critical surface, does not necessarily con- 
verge to a fixed point but to a periodic or even a nonpe- 
riodic orbit (we shall see that these attracting orbits are 
those of the Arnold's c at m ap) . Already if we start with 
the even perturbation (6.1) with M,P < 0, we obtain a 
cycle of period three as it has also been encountered in 
area- preserving maps [ p8[ . This can be understood by 
the symmetries of the transformation |]l9| , f29"|| . Starting 
with the two main resonances 1/1,1/2, the most general 
perturbation can be written as 

fo(<f) = - M e COS (l^i ■ (f) - P e COS (V> 2 • f) 

—M a sin (u\ ■ (f) — P a sin (i>2 • y) , (6-2) 
= -Mcos(i/ x • (ip + 9)) - Pcos{i/ 2 ■ {<p + 0)) , 



where M = ^Mf+Mj, P = ^P*+P*, 
and 

9 =(- arctan(M Q /M e ), arctan(M Q /M e ) - arctan(P Q /P e )) . 
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The question is to analyze the effect of a shift of the 
angles on the transformation. We define a shift operator 
by 



Tg : ip h-> if + 9. 



(6.3) 



The KAM transformation co mmutes wi th T g , as can be 
easily verified from Eq. or Eqs. (0-(fh7|). The 

action of Tg on the shift of the resonances N 2 is described 
by the following intertwining relation: 



TZ O Tg = T N 2g O TZ, 



(6.4) 



where TZ denotes the KAM-RG transformation. Ap- 
plying this relation to the fixed point H^(ip), we have 
KH^v + Q) = H*(ip + N 2 8). The following map (which 
is Arnold's cat map P0| ) 



i ► N 2 mod 2tt, 



(6.5) 



gives the nature of the orbit to which the transfo rma tion 
converges starting with the initial perturbation ( |6.2| ) on 
the critical surface. For instance, if we s tart with 6 — Ttv\ 
which corresponds to the perturbation (|6.l[) with M , P < 
0, the transformation converges to a cycle of period three 
because N 6 6 = 6 mod 2ir. More precisely this 3-cycle is 
the periodic sequence 



{H*((p + 7tvi + tt(vi - u 2 )),H*(ip + irv 2 )}. 



(6.6) 



For any orbit of the map (5.5), there is a fixed set that 
plays the same role as the nontrivial fixed point for the 
even perturbation (6.1). These invariant sets belong to 
the same universality class as the fixed point, and in par- 
ticular, they have the same critical exponents. 
We can define a modified renormalization transformation 
such that the KAM-RG transformation converges to the 
fixed point for all initial perturbations of the form 
( |6.2D on the critical surface, by adding an initial shift of 
the angles 71 g. 
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FIG. 1. Critical coupling e c (L) as a function of L, the 
size of the cell Cl containing (2L + l) 2 Fourier coefficients. 
The upper curve corresponds to the KAM transformation, 
and the lower one to the KAM-RG transformation. 



VII. CONCLUSION 



The results show that the KAM-RG transformation is 
a powerful tool to describe the breakup of invariant tori. 
In particular, the transformation describes with high ac- 
curacy the critical surface which is the stable manifold 
of a nontrivial fixed point (or more generally, of a non- 
trivial fixed set related to this nontrivial fixed point by 
symmetries). 

We have implemented the KAM-RG transformation for 
the torus with frequency vector ojq = (1/7,-1). The 
extension to other frequencies that are quadratic irra- 
tionals is relatively clear. The case of a general irra- 
tional frequency will involve qualitatively new features. 
The KAM-RG transformation we described was for sys- 
tems with two degrees of freedom, but the extension to 
three |H| or higher dimensional systems jl9| should be 
accessible. 
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FIG. 2. Critical coupling e c (L) of the KAM-RG 
transformation as a function of L, the size of the 
cell Cl containing (2L + l) 2 Fourier coefficients. 




FIG. 4. Weight of the Fourier coefficients of 
White: < 10" 10 , grey levels: [1CT 10 , 1(T 7 ], 
[10- 7 ,10- B ],[10- 6 ,10- 3 ], black: [lO -3 , 10" 2 ]. 




FIG. 3. Weight of the Fourier coefficients of /*: 
White: < 10- 10 , grey levels: [10- 10 ,10- 7 ], 
[10- 7 ,10- B ],[10- 6 ,10- 3 ], black: [10~ 3 , 10" 2 ]. 



FIG. 5. Weight of the Fourier coefficients of 
White: < 10- 10 , grey levels: [10" 10 , 10" 7 ], 
[10- 7 ,10- 6 ],[10- 5 ,10- 3 ], black: [lO -3 , 10" 2 ]. 
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